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Abstract. We present a numerical extension to the analyti- 
cal propagation model introduced in Hein and Spanier (2008 ) 
to describe the leptonic population in the galactic disc. The 
model is used to derive a possible identification of the com- 
ponents that contribute to the leptonic cosmic ray spectrum, 
as measured by PAMELA, Fermi and HESS, with an em- 
phasis on secondary e + — e~ production in collisions of cos- 
mic ray particles with ambient interstellar medium (ISM). 
We find that besides secondaries, an additional source sym- 
metric in e + and e~ production is needed to explain both the 
PAMELA anomaly and the Fermi bump, assuming a power- 
law primary electron spectrum. Our model also allows us to 
derive constraints for some properties of the ISM. 



1 Introduction 

Recently the leptonic component of the cosmic ray spectrum 
has gained new attention. New observations from ATIC, 
PAMELA and Fermi show a deviation from a power-law 
in the form of an excess in both the electron and positron 



. — spectra. Annihilating dark matter ( Allahverdi et al 



and nearby pulsars (Biisching et al. 2008 Grimani 



2009) 



2009 



Blasi and Amato 201 1|> (among other hypotheses) have been 



proposed as possible sources of the excess leptons. Regard- 
less of the source, a propagation model is needed to connect 
the energy spectrum measured on Earth with the injection 
spectra. 

We present our numerical cosmic ray transport model in 
application to the high energy electron transport in the ISM. 
Spatial and momentum diffusion, particle escape, accelera- 
tion via Fermi I and continuous energy losses were taken 
into account and their effects on the steady-state energy spec- 
trum analyzed. In solving the transport equation we em- 



ployed quasi-linear transport theory, the diffusion approxi- 
mation and a separation of the spatial and momentum prob- 
lem to obtain the leaky-box-equation, which was then solved 
numerically. The spatial problem was solved analytically in 
cylindrical and prolate spheroidal coordinates. 

The transport model was employed to calculate the spec- 
trum of secondary electrons in our galaxy. We assume the 
leptons from pioArXivn decay to be the dominating com- 
ponent of the secondary spectrum. The lepton spectrum 
from pion collisions of highly relativistic cosmic ray pro- 
tons with thermal protons in the ISM was calculated using 
the parametrization of pion production in p-p-collisions from 
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Kelner et al. (2006) and used as the injection spectrum for the 
transport model. With realistic simulation parameters the re- 
sulting positron flux in the vicinity of the solar system lies re- 
markably close to the low-energy PAMELA data. Assuming 
a generic power-law primary e~ spectrum and an additional 
e + — er source allows us to fit the datasets from PAMELA, 
Fermi and HESS very nicely and to put constraints on sev- 
eral transport model parameters and the corresponding ISM 
properties. 



2 The data 

The leptonic cosmic ray spectrum gained new attention when 
the results of the balloon experiment ATIC were published 
( |Chang et al.| |2008| l. The authors claimed an anomalous 
excess in the spectrum just below 1 TeV. This started wild 
speculations about the nature of the source for these leptons, 
which should be in the relative vicinity of the solar system. 

Those speculations were fueled further by new measure- 
ments from the PAMELA satellite published in |Adriani et| 
|aL] p009 ) claiming a rise in the positron spectrum above 10 
GeV. The lepton spectrum from the Fermi satellite (Abdo et 
al. 2009 1 could not confirm the ATIC-peak quantitatively but 
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showed a deviation from a straight power-law in the same en- 
ergy region. 

The HESS atmospheric Cerenkov telescope extended the 
spectrum measurement beyond 1 TeV showing a cut-off or a 
steeper power-law but could also not confirm the ATIC ex- 
cess ( |Aharonian et al.| |2009). 

Recently the Fermi lepton spectrum was extended to lower 
energies (7GeV) in Ackermann et al. (2010 1. Also the 



In this work we used the analytical solution for cylindri- 



PAMELA results were updated with new experimental data 



and different evaluation techniques in Adriani et al. (20101. 



To interpret these measurements, it is necessary to con- 
nect the injection spectra at the (possible) sources with the 
observations in the solar system. To do that we need a trans- 
port model that includes all the relevant processes in the ISM 
(spatial and momentum diffusion, cooling, particle escape, 
etc.). 



3 Transport model 

The derivation of the transport equation follows closely 
Lerche and Schlickeiser ( 1988 ) and Hein and Spanier (2008 ). 
Using Quasi Linear Theory and the diffusion approxima- 
tion the Vlasov equation is transformed into the steady-state 
diffusion-convection equation for the phase space density 
f(x.,p), which can be written in the form 



C x f+£ p f + S(x,p) = 
with the spatial operator 
£ x (x,p) = V[k(x,p)V] 
and the momentum operator 

\ d ( d 

£ p (x,p) = -— i [ p 2 D(x,p) — -p 2 p(x,p) 



r(x,p) ' 



(1) 



(2) 



(3) 



This equation describes the diffusion (coefficients k and D in 
space and momentum, respectively), acceleration and cool- 
ing (rates p) and escape (timescale r) of an injected source 
distribution S. 

Following |Lerche and Schlickeiser (1988) the equation 
can be solved by separating the operators and the source term 
in their spatial and momentum dependencies. Then the re- 
sulting steady state particle distribution can be written as an 
infinite sum: 



f(x,r) = YA i (x)R i (p) 



(4) 



The spatial coefficients Ai (x) can be acquired by solving the 
diffusion equation in the appropriate geometry 



8_ 

On 



A i exp(-wf-u) = V(kVAj6xp(-w i 2 'u)), 



(5) 



where wf take the role of inverse escape times for the mo- 
mentum modes and u is a convolution variable. 



cal geometries as provided by Wang and Schlickeiser ( 1987 1 
since this offers the best reflection of the symmetries of a 
spiral galaxy. Analytical solutions in spherical and prolate 



spheroidal coordinates are also available ( Schlickeiser et al. 



1987| Hein and Spanier 2008 1 and can be used to describe 



particle transport in elliptical galaxies and galaxy clusters. 

The momentum modes R t (p) are solutions of the ODE, 
commonly referred to as the leaky-box equation: 



l_d_ 
p 2 dp 



a 2 p 4 dR t aip 3 



k(p) dp n(p) 



R t +pRi 



■(^K(p)R i = -Q(p)(<S) 



The parameters a\ and a 2 represent the absolute scales for 
Fermi I and II processes and depend mainly on the Alfven 
and shock speed, respectively. As [Lerche and Schlickeiser| 
(1988) have shown, an analytical solution for the equation 
can be found, if the cooling rate p scales linearly with the 
particle energy, which is a good assumption for high en- 
ergy protons and heavier nuclei. Electrons above a few GeV, 
however, are cooled primarily by synchrotron radiation and 
IC scattering, both scaling quadratically in particle energy. 
Therefore the leaky-box equation was solved using numeric 
relaxation methods. 

To summarize, the model describes a particle distribu- 
tion that is injected by a homogeneous population of time- 
independent sources in the galactic disc (injection region), 
diffuses throughout the galactic disc and halo (confinement 
region), cooling by bremsstrahlung, adiabatic deceleration, 
synchrotron radiation and inverse Compton scattering, and 
leaves the confinement region at some point. The matter, 
photon and magnetic field distributions are assumed to be 
homogeneous over the whole confinement region. 

A model with a uniform source distribution is not a good 
choice to simulate high energy leptons from discrete sources 
(SNR, PWN. . . ), since the short energy loss timescales make 
the source distribution critical for the resulting spectra, as 
discussed in Blasi and Amato (2011). Conversely, produc- 
tion of secondaries from p-p collisions in the galactic halo 
can be approximated very well as a homogeneous source, 
since there are no dramatic variations in the CR and ISM 
densities. 



A similar analysis by Moskalenko and Strong ( 1998 ) used 
an analogue model as in this paper, but obviously did not 
include the new data above 1 GeV. Stephens ( 200 1 ) l has basi- 
cally the same drawbacks, in addition this work has a more 
sophisticated secondary production model. Compared with 
the recent works in this area like Blas i and Amato| ( |2011[ ), 
we concentrate more on the low energy end of the current 
PAMELA data and the constraints it can put on propagation 
parameters if interpreted as secondaries. 
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4 Secondary leptons 

Fortunately there is a lepton production process that can be 
approximated with a homogeneous spatial distribution: the 
secondary leptons from the decay of products of proton col- 
lisions. To get a galaxy-averaged flux of secondaries we 
consider collisions of high energy cosmic ray protons with 
ambient thermal matter, assuming that both populations are 
homogeneously distributed with average densities. 

The pion production cross-section and the electron source 
function for a single p-p collision used here are parametriza- 



ner et al. 



tions of results f rom the SIBYLL event generator from Kel 
(2006 1. Folding the source function F e ( jjf-,E, 
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dE 



with the cosmic ray spectrum J p {E p ) and the inelastic colli- 
sion cross-section ai ne \(E p ) yields the electron flux: 

cr incl (E p )J p (E p )F e (^,E e j -^ (7) 

Since the source function is sharply peaked, the CR spectrum 
can be cut off exponentially at the first knee (s=s 1000 TeV) 
with negligible deviations to the lepton spectrum below 10 
TeV: 

-2.677 



J P {E p ) = 0.252 
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This parametrisation of the CR spectrum was obtained by 

fitting the CR data between 10 10 and 10 14 eV in 

( fT997] l. 

The average thermal matter density uh is treated as a free 
parameter. The resulting electron flux with a slightly harder 
power-law (spectral index 2.62) is used as the injection func- 
tion in the transport model. Since the parametrizations in 



Kelner et al. (2006) only consider p-p interactions, the ac- 
tual flux of secondaries is underestimated by the contribution 
from heavier nuclei. 



5 Parameter constraints 

The principal shape of the resulting steady-state lepton spec- 
trum is shown in Fig. [T] The cooling, diffusion and escape 
processes introduce breaks in the injected power-law, so that 
multiple energy intervals with different spectral indices are 
formed, each of them defined by a different process. At 
the first break the injected power-law is steepened by 1/3 
as the diffusion/escape time scale becomes shorter than the 
time scale of bremsstrahlung and adiabatic cooling (linear in 
energy). The second break appears when the time scale for 
synchrotron and IC losses becomes shorter that the diffusion 
time scale. These quadratic losses steepen the spectrum by 
1 . The absolute flux value is determined mainly by the ISM 
density, since the injected secondary spectrum has a linear 
nn dependence (Eq. |7|), and to a lesser extent by escape 
and cooling losses, which transport particles out of the simu- 
lation region. 
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Fig. 1. Parameter studies for the secondary positron spectrum 
compared to the PAMELA e + fraction multiplied by a fit of the 
Fermi e + + e~ data. Combined systematic errors of PAMELA 
and Fermi are shown by the grey area. The thick black line with 
uh = 1.5- 10 3 m -3 , B = 3/iG and t — t CBC yields the best fit of the 
positron spectrum, if combined with an extra high energy source as 
shown in Fig. [2] The thin and dashed black lines show the slope 
change with changing escape losses. 



So the model has three major parameters: the linear energy 
loss time scale that reflects the thermal gas density nn, the 
quadratic loss time scale that incorporates the energy densi- 
ties of the EM field (composed of a constant CMB density 
of f=a 0.3 eV/cm 3 and a variable magnetic field B) and the 
escape time scale t csc that corresponds to the size of the con- 
finement region. 

Measurements of the positron spectrum, in particular the 
newest results from PAMELA, allow us to constraint those 
parameters. A rather conservative statement can be made, if 
we assume that the secondaries have a negligible contribu- 
tion to the positron spectrum. As shown in Fig. IT] assuming 
an average ISM density of sa 10 3 m -3 , the confinement re- 
gion has to be less than the assumed 30 kpc or the diffusion 
process has to be more efficient, so that the particles leave 
the galaxy 10 times faster (black dotted line). A higher ISM 
density would require an even more efficient particle escape, 
so as not to overshoot the measured positron spectrum. Since 
the spectral shape wouldn't play any role, the magnetic field 
can be almost arbitrary. 

On the other hand, the PAMELA data can be fitted very 
well with a superposition of two power-laws with a spec- 
tral index of about 3.6 for the low energy part, correspond- 
ing remarkably well with the secondary spectrum steepened 
by synchrotron and IC losses. Assuming the low energy 
positrons to be secondaries puts much stricter constraints on 
the parameters, since the second break in the spectrum has 
to appear at w 1 GeV to match the data. The break energy is 
determined by the equilibrium of the diffusive/escape losses 
and the radiative losses. The thick black line in Fig. [T] is 
the best fit with n H = 1.5 ■ 10 3 m~ 3 , B = 3fiG and t = t esc . 
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Fig. 2. A fit of Fermi (Ackermann 



et al. 20101 and HESS (Aharonian et 



E [GeV] 10000 



al. 2009) lepton data (solid red line) 
and PAMELA ( |Adriani et al~| |20l0l > 
positron data (solid green line). Grey 
areas show the systematic errors. The 
components that contribute to the spec- 
tra are secondary positrons (dashed 
green), a generic power-law primary 
spectrum (dashed red) and an additional 
symmetric e + — e~ source (solid and 
dashed black). 



The radiative losses in the galactic halo are dominated by IC 
scattering on the CMB photons, so that lowering the mag- 
netic field has no effect on the break energy. A slight flat- 
tening at the low energy end of the spectrum is needed for 
the best fit of the new PAMELA data. This corresponds to 
a halo radius of 30 kpc and a spatial diffusion coefficient 
n(E = lGeV) = 4.5 • 10 24 m 2 s _1 , yielding an average con- 
finement time i osc = 2 ■ 10 17 s. The absolute flux fixes the 
third parameter, the ISM density, to 1.5- 10 3 m -3 . Taking 
into consideration that about 99% of the confinement region 
represent the galactic halo with an average gas density of 
about 10 3 m~ 3 (McKee and Ostriker 1977 i, the parameters 
seem very realistic. 



6 Lepton excess 

In Fig. [2] the latest Fermi lepton measurement, in which the 
spectrum has been extended down to 7 GeV, is shown. A 
generic fit of it was used to calculate a positron flux from 
the PAMELA positron fraction measurement. This data is 
used because it was taken simultaneously in the same solar 
cycle making interpretation easier. The HESS high-energy 
data was shifted down by 15% to better coincide with the 
Fermi data. This is still within the systematic error margins 
of HESS, not to mention the systematic error of Fermi being 
of the same order of magnitude. 

Once the secondary background has been subtracted from 
the PAMELA positron flux, the excess can be fitted very well 
with a simple power-law with a spectral index of 2.3 (Pig. [2] 
black dashed line). Remarkably, if we assume the source of 
these "extra" positrons to be symmetric in e + — e~, we can 
also fit the Fermi bump leaving a primary lepton component 



of the shape E~ 32 • exp(— E/1A6T&V) (red dashed line). 
The (shifted) high-energy HESS data provides a means to fix 
the cut-off energy of the "extra" leptons to about 920 GeV. 



7 Conclusions 

With the help of our semi-analytic transport model we were 
able to identify the low-energy PAMELA data as secondary 
positrons from collisions of CR protons with ISM gas us- 
ing very realistic values for the simulation parameters. An 
indirect implication from that is that lepton transport is dom- 
inated by radiative losses (synchrotron and IC) above 1 GeV. 
Assuming a generic primary lepton spectrum of the 
form E~ 3 - 2 -exp(-E/lA6 TeV), the PAMELA anomaly, the 
Fermi bump and the high energy HESS data can all be fitted 
by a single component, that is symmetric in e + — e~ with 
a power-law spectral index of 2.3 and a cut-off energy of 
920 GeV. This could either be a lepton population injected 
into the ISM by nearby source with E~ 2 - 3 , so that the propa- 
gation time scale is shorter than the radiative loss time scale, 
or a more distant source with a much harder spectrum 25 -1 - 3 
which is then steepened by radiative losses. A supernova 
remnant would be a good candidate in the first case, pulsars 
of different ages would suite both cases. 
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